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INCOMPLETE  INTEGRALS  OF  CYLINDRICAL  FUNCTIONS 


INTRODUCTION 

The  class  of  cylindrical  functions  C  includes  Bessel  functions  of  the  first  kind  J,  modified 
Bessel  functions  /,  Bessel  functions  of  the  second  kind  or  Neumann  functions  T  (or  A'),  Bessel  func¬ 
tions  of  imaginary  argument  or  MacDonald  functions  K,  and  Bessel  functions  of  the  third  kind  that 
include  Hankel  functions  of  the  first  and  second  kind,  Hil>  and  H{2). 

The  general  incomplete  Lipschitz-Hankel  integral  of  cylindrical  functions  C„(z )  is  defined  as  the 
function  of  two  complex  variables: 

Ce Ja,z):  -  J  (1) 

Here  the  symbol  e  denotes  the  presence  of  the  exponential  function  and  n,  p  may  be  complex. 
Analogously,  we  define  integrals  that  contain  the  functions  sin  (at)  and  cos  (at)  in  place  of  exp  (at): 


cs "(a,  z):  —  sin  (at)  t^C^dt 

(2) 

Cc  (a,  z):  =  jo  cos  (at)  ttLCv(t)dt . 

(3) 

To  assure  convergence  of  Ce  (a,  z)  and  Cc  (a ,  z),  it  is  necessary  that  Re  (p.  +  1)  >  |  Re  r  |  when 
C  -  K,  Y,  H(l),  H(2) ;  Re  (1  +  /t  +  p)  >  0  when  C  =  /,  J .  For  convergence  of  Cs  (a,  z) 
replace  by  n  +  1  in  the  latter  two  inequalities.  When  fi  =  p,  we  define,  for  example,  Ce  =  Ce 
where  for  convergence  Re  n  >  —1/2  for  all  C. 

Integrals  of  the  type  given  by  Eqs.  (1)  to  (3)  occur  very  often  in  applied  mathematics.  Agrest 
and  Maksimov  [1]  have  found  representations  for  Ce(a,  z),  ,  (u  z),  and  Cc  (a ,  z)  using  incomplete 

cylindrical  functions  which  are  themselves  expressed  as  integra  '  n  this  report  which  is  an  extension 
and  revision  of  previous  work  [2]  we  give  representations  for  Ce  (a,  z),  Cs  (a,  z),  and  Cc  (a,  z) 
using  only  the  Kampe'  de  Feriet  double  hypergeometric  functions  F’-j-ijol*,  >’  ]  and  ,  y). 

PRELIMINARY  RESULTS  AND  DEFINITIONS 

To  begin,  we  summarize  some  results  that  are  found  in  Ref.  3,  p.  85:  Let  a  and  b  be  arbitrary 
constants, 


F„(z):  =  a/,(z)  +  beiv'Kr(z) 
G„(z):  =  aJr(z)  +  bYv(z) 


0:  = 


fl:  H  =  F 
(0:  H  =  G 
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Then 


-—•'ft* 


J  ^H,(r)Jr  =  e  2  [Ox  +  v  -  l)zH„(z +  (20  -  l)azH,_  ,(z)sM  .,(<*)] ,  (4) 
where  the  Lommel  functions  „  Ox  ±  *  *  -1,  -3,  -5,  . . .)  are  given  by 

•V  „(2)  = 


T^+l 


Ox  -  V  +  1)0*  +  V  +  1) 


lr  2 


H  ~  v  +  3  ii  +  v  +  3  -z2 

2  2  ’  4 


(5) 


Now  defining 


*:  = 


r 

r  r 

1:  C  =  I,K 

1:  C  =  K 

—  1:  C  =  //,  J,  y  7;  = 

-1:  C  =  H,I,J,Y  w:  =  “ 

L  l 

1:  C  =  AT 
2/ir:  C  ~  H ,  Y 
0:  C  =  /,7 


we  may  deduce  from  Eqs.  (4)  and  (5)  the  result 


.2  +  1 


Jn^C'fCO  *  = 


H  -  v  +  1 
_  ,(z) 


*\(z)  ,F2 


l; 


.  M  -  t>  +  3  n  +  V  +  1  £z2 


’  4 


fi  +  v  +  \ 


lr  2 


l; 


.  M  -  y  -f  3  II  +  v  4  3  fc2 


(6) 


We  define  the  Kampe  de  Feriet  double  hypergeometric  functions  L,  M  and  Q  and  give  associ¬ 
ated  generating  relations  [4,  5]: 


L[ or,  0;  7,  5;  x,  yj: 


_  IT  0: 1 ;  1 
*  2:0;0 


M[a,  0;  7,  5;  *,  y]:  =  F  }:?$ 


- :  a;  0; 


7,  6:  - ;  - ; 


y 


\x  I  <  OO  ,  \y  I  <  00 


Q[a,  0,  7;  il,  v,  X;  x,  y]:  = 


0; 


x,  y 


7:  5;  - ; 

- :  &\  Y, 

L  M,  v:  X;  - ; 


x,  y 


\x  I  <  00  .  |y  |  <  l 


|  x  |  <  00  ,  |  y  |  <  00 


“  (oOm  Xm 

L[a,  0;  7,  6;  x,  y]  =  £  ,  —j  m  +  7>  m  +  5;  y] 

m  =0  (7)mWl  m- 


(7) 


2 


(8) 


M[a,  /3;  7,  6;  x,  or]  =  £ 


(<*)m  xm 


7  3^o[0.  1  ~  6  -  nr, /] 


=0  (7)«(«)*  »«! 

00  (a)-(fi)~  rm 

Q[a,  0,  y;  fi,  v ,  X;  y]  =  £  7  7  /  - j~  1^2 ™  »  +  *;  >’3 

m=  0  (#»)«(•')*  W*  m! 

It  is  easy  to  see  that  the  function  L  is  a  special  case  of  Q : 


For  brevity  we  define  the  parameter  lists 


A 1  (ji,  v):  = 

A2(ji,  v):  = 

B\in,  v ):  = 

B2(ji,  *>):  = 

D.Ox):  = 

D2(p)'-  = 

E i0*,  k):  = 

£20,  v):  = 


2 

2 

H  +  v  +  1 

H  —  v  +  1 

2 

2 

/t  +  v  +  2 

H  -  v  +  2 

2 

2 

n  +  V  +  2 

H  —  f  +  2 

2 

2 

y  +  i; 

1  .  3 

2  +«.y 

—  +  i; 

3  4-  3 

2  ^  2 

1  +  ft  +  V 

2  +  n  +  v 

2 

2 

1  +  n  +  p 

2  +  n  +  v 

2 

2 

,  l; 


,  l; 


,  l; 


,  I; 


]  =  £[a,  0\  y. 

&>x,  y]. 

H  +  v  +  1 

H  —  v  +  3 

2 

2 

n  +  V  +  3 

H  —  v  +  3 

2 

2 

H  +  v  +  2 

H  —  v  +  4 

2 

2 

H  +  v  +  4 

H  —  v  4-  4 

2 

2 

GiO*)'-  =  1  + 

1 

1  ♦  2  ’  ^ 

GzOO-'  =  1  + 

1;  y>  1 

2  +  n  +  v 

3  +  (i  +  »» 

2 

2 

3  +  n  +  v 

4  +  ft  +  v 

’  2 


,  1  +  * 


,l  +  v 


F\(n):  =  |  +  1;  1  +  M,  |  +  M  W,0t):  =  2  +  ji,  1,  1;  2  +  ft,  3, 


^O1):  =  y  +  F.  l;2  +  ft,  y  +  p  #2(**):  =  2  +  n,  1,  1;  3  +  p,  3, 


(9) 


3 


N>|C*  K»|c* 


REPRESENTATIONS  FOR  Ce  (a,  z).  Cs  (a,  z),  Cc  (a,  z) 

Substituting  the  Maclaurin  series  for  exp  (at)  in  Eq.  (1)  and  splitting  into  even  and  odd  terms 
we  obtain  on  integrating  term  by  term 

c'- z)  ' i« " +  +  .£  TTT^  J^'1  *  (10) 

Then  using  Eq.  (6)  and  the  generating  relation  Eq.  (9)  we  obtain  after  a  tedious  but  straightforward 
computation  the  result 


Ce  (a,z)  =  zl+>C,(z) 


r  ■  c 

.  a2z2  £z2 

l  ^  iO 

R  .  «2Z2 

id' 

1 

Im-J'  +  I  Q 

+  1  0  V 

H  —  v  +  2 

[B"  4 

’  4 

J 

+  VZ2  +  >Cv^(z) 


r 

r  .  «V  fe2 

|^  +  p  +  l)(/1-*'  +  l)  * 

[Al-  4  '  4  J 

az 


(/!  +  »<  +  2 )(ji  -  v  +  2) 


r\  .  ah 2  iz2 

a  ’  —r~ 
4  4 


(id 


Since 


we  may  write 


Q.  S°  >  Z)  ~  ^7  j  Ce.  MQ'  Z )  “  Ce,  S~ia'  Z ) 


Q,  (a  ,  z)  =  —  -I  Ce>  (ia  ,  z)  +  Ce  (- ia  ,  z) 


Q  (a,z)  = 


az 


2+m 


H  —  p  +  2 


Cr(z)fi 


„  .  -a2z2  £z2 

r~,_r 


+7^c--'we 


r  iZ2V  fei1 
®2’  4  ’  4 


(12) 


Cc  (a.z)  = 


,i+n 


H  —  v  +  1 


C,(z)fi 


.  -a2z2  {z2 

:  -  r- 


_nz 


(l  +  v  +  1 


C,-i(z)G 


,  ~a2z2  id 

^  4  ’4 


(13) 


4 


For  n  =  v,  Eqs.  (11)  to  (13)  reduce  to 


Ce  (a,  z)  =  z 


=  T1  + 


'C/z) 


n.£^ 1  fel 

4  ’  4 


+ 

2 


a2z2  £z^ 


i0*.  m); 


+  iz2  +  'C,_l(z)  j— Z. 


02; 


flV 

4  ’  4 


+  —s_e 

4(1  +  m) 


z? 20*’  #*); 


4  ’  4 

a2Z2  £z2 
4  ’  4 


Q  (a,  z)  =  —  az 


1  2  +  u 

=  -  /I7*  ^ 


CJz)Q 


— a2z2  fz2 


vz 


2(1  +  M) 


C„_i(z)G 


— a2z2  £z2 

«**.  <*>;  ~x~  •  4 


Cc(a,z)  =  z1+MCM(z)L 


z>,; 


2  2 


fei 

4 


VZ 


1  +  2/i 


C,_,(z)L 


z>2; 


-a2z2 


fei 

’  4 


Defining /  +  :  =  J,J~:  =  /,  it  is  interesting  to  note  that  we  may  also  write  [7] 

az 


J*  (a  ,z)  — 


z  \  +ii  +  »eaz 

2\\  +  ^  +  p)r(i  +  v ) 


Eu 


TZ2  a2z 2 
4  ’  4 


2  +  fi-+  v 


Q 


(14) 


(15) 


(16) 


E 2; 


tz2  a2z2 


4  ’  4 


•/,*  (a,  Z)  = 


z(z2/2YeK 

L .  ^  -v 1 

oz 

.  =Fz2  a2z2 

1 

(i  +2foro  +#o  \ 

F"  4  ’  4 

2(1  +m) 

2’  4  ’  4 

J 

(17) 

(18) 


Here  the  Bessel  functions  J  *  do  not  appear. 

We  remark  that  Eqs.  (11)  to  (18)  are  valid  for  ja|  <  oo,  |z|  <  oo. 

OTHER  REPRESENTATIONS  FOR  Ce  (a,  z),  C5  (a,  z),  Cc  (a,  z) 

In  Eq.  (4)  sM  „  may  be  replaced  by  the  Lommel  function  „  where,  when  either  of  the 
numbers  n  ±  v  is  an  odd  positive  integer, 


V  r(Z)  =  Z' 


2/E  0 


1, 


v  —  n  +  1 


-  1 


-4 

,2 


(19) 
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Defining 


1 

(-1V2/T 

0 


C  =  K 
C  =  H ,  Y 
C  =  1,J 


and  replacing  „  in  Eq.  (4)  by  v  given  by  Eq.  (19)  we  obtain  for  nonnegative  integers  n  the 
result 

J  ‘  f  +  2n  +  1  Cv(t)dt  =w„2>  +  2n  T(n  +  1)  T {v  +  n  +  1) 

-  fe”  +  2n  (  2(v  +  n)  C„(z)  jF ol  1  >  -n,  1  -  v  -  n\— ;  4£/z2] 


+  ijzC„  _  ,(z)  jFqII,  -n,  -v  -  n ;  — ;  4£/z2]  ). 


(20) 


Using  Eqs.  (6)  and  (20)  together  with  Eq.  (10)  and  noting  the  generating  relations  Eqs.  (7)  and 
(8)  we  arrive  at 


Ce  (a,  z)  ~  a>  2M  T(1  +  n)  a  2F\ 


1,  l  +  y;  &2 


(21) 


+  z*C  (z)  \z  L 


D\  (m); 


fe2 


-  2 ZapM 


G,(4);~^,^2 


+  72  1  +  '‘C  l(2) 


1  +  2/i 


n  a2z2  fc2 
— .  -j- 


-  &  M 


,2,2 


G2(m);  ^ .  fr2 

4 


|  a  |  <  1 ,  Re  n  >  -  1  /2 . 

Since  this  computation  is  straightforward,  it  is  omitted. 

In  order  to  obtain  formulas  for  the  C,o(a ,  z),  i.e.  for  a  —  0,  from  Eq.  (21)  we  may  readily 
show  that 


ax 


lim  5  M  [a,  0;  7,  6;  *,y]  =  —  M  [a  +  1,  j8;  7  +  1,  2;  x,  y] 
6-0  7 


from  which  it  follows  that 


lim  2/i  M 
v-  -  0 


G,00;  ~~ ,  fa2 


y  a2z2  W 


«  „2„2 

2,  1;  |,2;  £a2 

2  4 


6 


We  also  note  that  when  n  =  0  in  Eq.  (21), 


so  that  also 


2^1 


sin  la 
1  —  a2 


sinh^a 
a'J  1  +  a 2 


The  result  obtained  from  Eq.  (21)  for  Ke  {a ,  z)  is  given  in  [4]. 

On  setting  n  =  v  in  the  equations  for  Cs  and  Cc  given  in  the  previous  section  in  terms  of 
Ce  we  obtain  on  using  Eq.  (21) 


Cs(a ,  z) 


co  a  2^  T(1  +  H)  2^1 


1. 1  +  y;  -la2 


(22) 


-  Zaz"  i  2fiCJz )  M 


2  2 

^  .  —a  z  fc_3 
— : — .  -la 


+  VZCp  ~  i(z)  M 


r  ~°2z2  >  2 

C2;  — t— - 


|a|  <1,  Re  jc  >  - 1 . 

For  Cc_(a ,  z)  we  obtain  again  the  result  given  by  Eq.  (16). 

It  is  shown  in  Ref.  5  that 


M[a,  1;  7,  5;  tx,  r]  =  1  +  qFJ-;  6;  *]{2F,|a,  1;  7;  r]  -  1) 


V 


c*tr 


278(6  +  1) 


Q[  1  +  a,  1,  1;  2  +  8,  3,  1  +  7;  tx,  jc J . 


(23) 


This  equation  provides  the  corollary  that  M[a,  1;  7,  8;  tr,  r]  converges  on  the  unit  circle  1 1  |  =  1  if 
and  only  if  2Fi[a,  1;  7;  r]  does.  We  then  deduce  that  Eqs.  (21)  and  (22)  are  conditionally  conver¬ 
gent  on  \a  |  =  1,  ±  fa2  *  1  provided  that  (Re/c|  <  1/2. 


We  may,  however,  use  Eq.  (23)  to  better  advantage.  We  find 


M 


fjZ.2 

Gl(");  4  >  &2 


=  1  + 


V'"* 


r(p)j~il{z)\2Fl 


1  +  v,  1;  y;  la2 


-  1 


|a2(z2/4)2 

3p 


2  2 
a  z 


lL2 

4 


7 


M 


r  t  \  °2z 2  t  2 

G2(&,  —r~ .  & 


L 


=  1  + 


ru  +  v)  j;hz)  \  2fx 


1  +  v,  1;  ta2 


«fl2fe2/4)2 


H2(v)\ 


,2,2 


_  id 

4  ’  4 


3(2  +  p) 

where  we  agree  that  J  * 1  =  J  * . 

Using  these  equations  together  with  Eq.  (2i)  and  noting  the  Wronskian  relations 

+  i  (z)/„(z)  +  K,(z)Iv  +  i(z)  =  1/z 


-  1 


J „  +  i(z)T„(z)  -  •/„(z)>'„  +  ](z)  =  2/irz 


we  deduce  after  some  computation 


Ce  (a  ,  z )  =  0)2'  TO  +  p)a 


(24) 


+  z,'C,(z)i  z  L 


„2,2  t,2 

W'^T-^r 


+  ■~T~~  0 

24 


id 

4 


—  2£a  p 


+  i,Z 1  +  rCr  _  ,(z) 


1  +  2p 


o2(»'); 


a2z 2  id 

4  ’  4 


a3z4 


48(2  +  p) 


H2(v)\ 


2  2 
a  c 


id 

4 


-  £a 


And  in  the  same  way  we  obtained  Eq.  (22)  we  find  from  this  result  that 

Cs  (a ,  z)  =  u  2”  T(1  +  v)a 


\ 


2  2  y  2 

-a  z  £z 


HAv)\ 


+  2iv 


TfdZ 


1  +  "C 


i>  -  1 


(Z) 


r  a 


[48(2  +  p) 


*»  T 


H2(v)-  ^ 


|r 

4 


+  i 


(25) 


INCOMPLETE  WEBER  INTEGRALS 

The  incomplete  Weber  integrals  Cw  (x,  z)  are  defined  for  Re  v  >  -1,  1 2X  -  arg  x  |  <  ir/2, 
z*:  =  a>  exp  /X  as  follows: 


J».Cr,  z):  =  J  w(x,  p;  t)  J ,{t)dt 


8 


Iw  (x,  z ):  =  Jo  W(x,  V,  t)  I „(t )dt 


Yh.  (x,  z):  =  j  "  w(x,  v,  t)  Y „(t )dt 

Kw(x,  z):  =  f  “  w(x,  v\  t)  Kv{t )dt 
'  Jz 


where  the  subscript  w  indicates  the  presence  of  the  kernel  defined  by 

w[x,  v,  t ):  =  (/ /2x)v  +  1  exp(;t  -  t2/Ax). 

It  is  shown  in  [1,  p.  136]  that  CWo(x ,  z)  may  be  given  in  terms  of  Cfo(a,  z).  Further,  Cw  (x ,  z) 
is  completely  determined  by  CK.o( x,  z )  when  the  index  v  is  an  integer  via  a  recursion  relation 
[1,  p.  123].  Hence,  CH; (x ,  z)  may  be  given  in  terms  of  Ce  (a  z)  or  equivalently  in  terms  of  the 
Kampe'  de  Feriet  functions  L ,  M  or  L ,  Q .  The  interested  reader  can  easily  produce  formulas  for 
Cw  (x,z)  in  terms  of  these  Kampe'  de  Fe'riet  functions.  We  remark  only  that  the  definition  of 
Ceo(a ,  z)  given  in  [1,  p.  135]  uses  exp(-ar)  in  place  of  exp(ar)  which  is  used  in  Eq.  (1). 

REDUCTION  FORMULAS  FOR  L  AND  Q 

Many  special  cases  of  Eqs.  (11)  to  (18),  (24)  and  (25)  may  be  obtained  in  one  form  or  another, 
provided  we  know  a  reduction  formula  for  either  L  or  Q.  We  summarize  a  few  known  relevant 
reduction  formulas  [4-8]: 


L[a,  0;  7,  6;  z,  z]  =  ,F2[q:  +  /3;  7,  6;  z] 


\ 


sinh  z 
z 


L 


2/x  sinh  z  cosh  z 
1  +  2/t  z  l+2/i 


Q 


B2(p-,  #t); 


1  +  /i  _4_ 
l+2/i  z 2 


cosh  z  - 


T(1  +  fi)I  „(z ) 


Q 


/1); 


2  ± 
l+2/i  z 


+  sinh  z 


ra  +  /i  )/„  _ 


Other  properties  and  reduction  formulas  for  L  and  Q  re  found  in  Refs.  4-8.  In  particular,  reduction 
formulas  for  Q]Ax{ji,  v)\  x ,  x],  Q[A2(ji,  v)\  x ,  x],  Q[B x(ji,  v)\  x ,  x]  and  Q[B2(ji,  i>);  x ,  x]  are 
derived  in  Ref.  8. 
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APPLICATIONS 


Of  interest  in  applications  are  the  Junctions  Jto(a,z),  le(a ,  z),  Ye(a  ,  z),  and  Kf  (a  ,  z). 
Je(a,  z)  and  Yto(a ,  z)  occur  in  problems  in  the  theory  of  diffraction  in  optical  apparatus  [1,  p.  227], 
The  function  Ie^(a ,  z)  plays  an  important  role  in  the  study  of  oscillating  wings  in  supersonic  flow  and 
arises  in  the  study  of  resonant  absorption  in  media  with  finite  dimensions  [1,  p.  195],  Kei,a,z) 
occurs  when  the  statistical  distribution  of  the  maxima  of  a  random  function  is  applied  to  the  amplitude 
of  a  sine  wave  in  order  to  calculate  the  distribution  of  its  ordinate.  This  latter  distribution  is  of 
interest  in  the  study  of  the  scattered  coherent  reflected  field  from  rough  surfaces  [9], 

Incomplete  Weber  inte^.als  of  order  zero  occur  for  example  in  the  theory  of  exchange 
processes,  the  analysis  of  drift  carriers  in  semiconductors,  and  high  temperature  plasma  physics.  See 
Ref.  1  for  further  details  and  other  applications. 

Since  the  functions  Ce(a ,  z )  are  of  some  importance,  by  using  Eq.  (14)  and  defining 


L^x,  y)  : 


L[D ,(0);  x,  y] 


L0(x,y):  =  L 


1  1  2. 

2  ’  l'  2  ’  2'X'y 


L[D2(oy,  x,  v] 


Qx(x,y):  =  Q 


1,  1,  I;  1,  2,  ~;x,y 


=  Q[Bx(0,  0);  x,  y] 


Qo(x,  >')■■  =  Q 


1,  1,  1;  2,  2,  — ;  x,  y 


=  Q[B2(0,0)-x,y] 


we  obtain 


Ke{,o,  z)  =  zK 0(z)  i  L j 


r  r 

r 

a2z2  z2 

L 4  ’  4  j 

+  -e. 

2  1 

a2z2  z2 

4  ’  4 

] 

+  z2Kx(z)  i  Lq 


r 

r  -\ 

L 

a2z2  z2 

Lo 

4  '  4 

^  J 

+  fQo 


a2z2  z2 


4  ’  4 


Yto(a,  z)  =  zY0(z) 

+  z2T,(z) 
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J'Sa'  z)  =  z^o(z) 

+  Z2J\{z) 

Ito{a,z)  =  z/0(z)  ■ 

-  z2Ii(z) 


The  equations  for  H}1)  and  H®}  are  the  same  as  those  for  Y(o  or  Je<>  with  Y  or  J  replaced  by  Hw  or 
H(2).  Further,  from  Eq.  (18)  we  have 


J'Sa '  z)  =  zeK 


,  I.  i.  shL  -zzi 

2  2  4  4 


L 


i  !•  2  2' 

2  2  4 


a2z2  -z2 


4o(a,  z)  =  ze“ 


1.  2  i.  fl2z2  ii 

2  ’  2  ’  ’  4  ’4 


Here  we  have  used  the  properties  of  L  that 


L[a,  0;  7,  6;  x,  y]  =  L[a,  0;  6,  y;  x,  y]  =  L[0,  a;  7,  6;  y,  x]. 


The  latter  results  for  Ce  (a,  z)  should  prove  useful  in  numerical  computation  of  these  functions.  We 
remark  that  other  formulas  for  Ce(,a,  z)  may  be  obtained  from  Eqs.  (21)  and  (24). 

SUMMARY 

Representations  for  incomplete  Lipschitz-Hankel  integrals  of  cylindrical  functions  using  only  the 
Kampe  de  Feriet  functions  in  two  variables  L,  M ,  Q  are  given.  In  addition,  several  known  relevant 
reduction  formulas  for  these  functions  are  provided.  Further,  it  is  indicated  that  incomplete  Weber 
integrals  of  integer  order  may  be  expressed  in  terms  of  incomplete  Lipschitz-Hankel  integrals  of  order 
zero.  Hence  incomplete  Weber  integrals  of  integer  order  may  also  be  represented  by  the  functions 
L,M,Q. 
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